How to measure squeezing and entanglement of Gaussian states without homodyning 
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We propose a scheme for measuring the squeezing, purity, and entanglement of Gaussian states 
of light that does not require homodyne detection. The suggested setup only needs beam splitters 
and single-photon detectors. Two-mode entanglement can be detected from coincidences between 
photodetectors placed on the two beams. 
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The recent rapid development of quantum information 
theory has largely stimulated research on nonclassical 
states of light, with the main focus on the generation of 
entangled states of light that are required for tasks such 
as quantum teleportation, dense coding, or certain types 
of quantum key distribution protocols. A particularly 
promising approach consists in processing quantum in- 
formation with continuous variables [l| , where the quan- 
tum information is encoded into two conjugate quadra- 
tures of the quantized mode of the optical field. The 
main advantage of this approach is that many protocols 
can be implemented by processing squeezed light into 
linear optical interferometers followed by measurements 
with highly efficient photodiodes Such experiments 
can be described in terms of Gaussian states which thus 
play a central role in continuous-variable quantum in- 
formation processing. In particular, squeezed Gaussian 
states provide the necessary source of entanglement. The 
squeezing is usually observed with the use of a balanced 
homodyne detector, where the signal beam is combined 
with a strong local oscillator (LO) providing a phase ref- 
erence . The observed quadrature fluctuations depend 
on the relative phase between the LO and the signal. The 
maximal squeezing of the signal then corresponds to the 
minimal observed quadrature variance. 

Given that quadrature squeezing is inherently a phase- 
sensitive phenomenon, one would expect that it may not 
be possible to determine the squeezing properties without 
an external phase reference (LO). In this paper, we show 
that, surprisingly, a phase- insensitive device is sufficient 
provided that we can a priori assume that the optical 
mode is in a Gaussian state. The setup we suggest con- 
sists in beam splitters with variable splitting ratios, phase 
shifters, and photodetectors with single-photon sensitiv- 
ity (e.g., avalanche photodiodes). It can also be extended 
to estimate the squeezing of multimode Gaussian states. 
In particular, a variation of our setup is capable of mea- 
suring the degree of entanglement of a two-mode Gaus- 
sian state, namely the logarithmic negativity [HQ. Be- 
sides the determination of squeezing and entanglement, 
our setup can also be used to measure the purity of Gaus- 
sian states . In addition, the detectors need not be per- 
fect and an efficiency rj < 1 can easily be compensated 



by proper data processing. 

Our scheme works for an arbitrary number of modes 
N and is economic with respect to N in the sense that 
the number of measured parameters is only linear in N 
while the full tomography of Gaussian states revealing 
the whole covariance matrix would require the measure- 
ment of oc N 2 parameters. In this context, it is re- 
lated to several recent proposals on how to directly mea- 
sure the purity, overlap, and entanglement of quantum 
states without full state reconstruction 0, S S 01 ■ It is 
also reminiscent of the photon-number distribution mea- 
surement scheme using a photodetector without single- 
photon resolution as proposed in Ref. [Icj . 

Preliminaries. Let us begin with introducing the 
necessary notation and definitions. Let r = 
(x%,pi, . . . , Xk,Pk, ■ ■ ■ , %n,Pn) be the vector of conjugate 
quadratures of N modes which satisfy the canonical com- 
mutation relations [cCj,pfc] = iSjk- The Gaussian state is 
fully described by the vector of mean values £j = (rj) 
and the covariance matrix 



(Ar 3 Ar k ) + (Ar k A rj 



(1) 



where Arw 



£j. The quantum state of the optical 



field can be fully characterized by a s-parametrized qua- 
sidistribution which provides phase-space representation 
of the quantum state. For our purposes, it is convenient 
to utilize the Husimi Q-function. The Q-function of an 
TV-mode Gaussian state is the Gaussian distribution 
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Q(r) 



v/det(7 + 7) 



exp[-(r-0 T (7 + ir 1 (r-0] 



(2) 

where I is the identity matrix. The squeezing properties 
do not depend on £ and are fully described by 7. The 
maximal observable squeezing, i.e. the minimal quadra- 
ture variance, is called the generalized squeezing variance 
A and can be determined as the minimal eigenvalue of the 
covariance matrix , 



A = min[eig(7)]. 



(3) 



The purity of the mixed state with density matrix p is de- 
fined as V — Tr[p 2 ]. For a Gaussian state with covariance 
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FIG. 1: Direct measurement of a Gaussian state with the 
use of a photodetector PD with single-photon sensitivity, pre- 
ceeded by a beam splitter BS of transmittance T. (a) Single- 
copy scheme for £ = 0. (b) Two-copy scheme for £ 7^ 0. 



matrix 7, one obtains [3 

V = [det( 7 )]- 1 / 2 . 



(4) 



Single-mode case. For the sake of simplicity, we first 
illustrate the procedure on single-mode Gaussian states 
(N = I). Consider the setup depicted in Fig. 1(a). The 
input mode impinges on a beam splitter BS with tunable 
transmittance T, and the output mode is measured by 
a photodetector PD with efficiency rj that is sensitive to 
single photons (no single-photon resolution is needed). 
We assume that this realistic detector can be modeled 
as a beam splitter with transmittance n followed by an 
ideal detector that performs a dichotomic measurement 
described by the POVM elements n = |0)(0| and ITi = 
1 — n . (In what follows, we assume that the detector 
is ideal and rj can be taken into account by substituting 
T — > T)T.) The probability of no-click of an ideal detector 
PD is given by P = Tr[pIT ] = (0|p|0) = 2ttQ(0), so that 
inserting r = into Eq. (JSJ yields 



Vdet( 7 ' + 7) 



exp[-r(7' + /)-Y] 



(5) 



where 7' and £' are, respectively, the covariance matrix 
and the displacement vector of the beam impinging on 
the photodetector. 

Suppose that we set the beam splitter transmittance to 
the value Tj . The covariance matrix 7' of the state after 
passing the beam splitter reads 7' = Tjj + (I — Tj)I. 
Similarly, the coherent signal is damped to £' = yTj£. 
On inserting 7' and £' into Eq. (JSJ, we obtain 



P, 



[det(7T,)] 1/; 



• cxp (-7X T 7t/£ 



(6) 



where jt- = Tj7 + (2 



Tj)I- 



We thus find that P,- 



depends on Tj and four parameters of the state: det 7, 
Tr( 7 ), 9i = C?722 + e 2 2 7n - 271266, and g 2 = £ 2 + & 
This immediately suggests that if we measure Pj for at 
least four different Tj's, we might be able to reconstruct 
the values of these four parameters by solving a system 
of nonlinear equations. However, numerical simulations 
reveal that the inversion of these highly nonlinear equa- 
tions typically leads to extremely large fluctuations of 



the estimated parameters even for a very large number 
of measurements for each setting Tj . 

Fortunately, in the important case where the displace- 
ment vector is zero (£ = 0), the scheme provides reliable 
and well-behaved estimates of these parameters since for- 
mula JfjJ) then simplifies to 

det 7 T . = T 2 det( 7 ) + T(2 - T)Tr( 7 ) + (2 - T) 2 = 4Pr a . 

(7) 

This results in a system of linear equations for det (7) and 
Tr(7). If measurements for two different transmittances 
T\ and T 2 are performed and the observed probabilities 
of no-click are P\ and P2, then the system of Eqs. (Q) 
can easily be solved and yields 



Tr( 7 ) = 



T2-T1 VTiPf 



To 



T 2 P? 



2 2 
2----, (8) 
J-i ±2 



det(7) 



T\—Ti \ TiPf 



2 - T 2 2-TA (2 - Ti)(2 - T 2 ) 



T{T 2 



Let us investigate what can be extracted from the knowl- 
edge of Tr(7) and det (7). As noted above, the squeezing 
properties of the Gaussian state, namely the generalized 
squeezing variance A, can be determined from the eigen- 
values of 7, cf. Eq. J2J. For a single-mode state, 7 is 
symmetric 2x2 matrix and its eigenvalues can be ex- 
pressed in terms of det (7) and Tr( 7 ), which are both 
determined by the present method. We find that 



Tr( 7 ) 



4 det (7) 



(10) 



Moreover, since our method provides an estimate of det 7, 
we can also determine the purity from Eq. 

If £ 7^ 0, our scheme is still usable provided that we 
can perform a collective measurement on two copies of 
the state, as depicted in Fig. 1(b). The two input modes 
Aj and A lr prepared in identical Gaussian state interfere 
on a balanced beam splitter. In the Heisenberg picture, 
the annihilation operators of the output modes A + and 
A_ are linear combinations of those of the input modes, 
a± = 2~ 1 / 2 (ar±aj/). The covariance matrix of mode A_ 
is equal to 7 [ljj but the coherent signal in A^ vanishes 
due to the destructive interference, £_ = 0. As shown in 
Fig. 1(b), the mode A^ is subsequently sent to a direct 
measurement setup identical to that shown in Fig. 1(a). 

Multimode case. We now extend this procedure to 
multimode Gaussian states. A reliable operation again 
requires two copies Aj and An. Essentially, we use 
in parallel N setups such as shown in Fig. 1(b). We 
combine each pair of modes Ai t k an d Azr,fe) with k = 
1,...,N, on a balanced beam splitter. Each "minus" 
mode j4_,fc, with £ = 0, is then sent to a beam splitter of 
transmittance Tj followed by a photodetector. We mea- 
sure the probability Pj = 2 N / (detj Tj ) 1/2 that none of 
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the N detectors clicks. For an iV-mode state, the deter- 
minant of 7t = Tj + (2 — T)I can be expanded as 

2N 

detftr) = £ T "( 2 - Tf N - n Ul) + (2 - T) 2N , (11) 

n=l 

where f n (j) is an homogeneous polynomial of nth or- 
der in the matrix elements of 7, e.g., /2jv(t) = det(7) 
and /i(7) = Tr(7). The probability Pj thus depends on 
Tj and the 2N parameters f n {j). If we measure Pj for 
2iV (or more) different transmittances Tj's, then we can 
determine the parameters /„ of the Gaussian state by 
solving a system of linear equations 

2N 

Y J T?{2-T 3 ) 2N - n f n { 1 ) = 2 2N Pr 2 -{2-T 3 ) 2N . (12) 

71=1 

Once we know f n {'~f)i we can determine the generalized 
squeezing variance A as the smallest root of the charac- 
teristic polynomial det (A/ — 7) = 0. It can be seen from 
Eq. (jl 1|) that the parameters fj (7) are the coefficients of 
this characteristic polynomial, and we have 

2N-1 

X 2N + A fc (-l) fe ./W- fe (7) = 0. (13) 

We can also determine the purity of the A^-mode Gaus- 
sian state from J2n(j) with the help of formula (@J. 

Entanglement detection. In the context of quantum 
information processing with continuous variables, the en- 
tanglement properties of Gaussian states deserve partic- 
ular attention. It has been shown that a two-mode Gaus- 
sian state is separable iff it has a positive partial trans- 
pose ^1 ^| . This property can easily be checked if one 
knows the covariance matrix 

1AB = ( V ^ ) (14) 

of the bipartite state, where ja and 7s are the covari- 
ance matrices of modes A and B, respectively, while gab 
captures the intermodal correlations. Moreover, analyt- 
ical formulas for several entanglement monotones that 
measure the entanglement of Gaussian states have been 
given in the literature 0, \\(\ . A particularly simple for- 
mula has been obtained for the logarithmic negativity 
Ej\f of an arbitrary Gaussian state. To calculate Ej\f, we 
must determine the symplectic spectrum of the covari- 
ance matrix of the partially transposed state p^s- As 
shown in Ref. , the symplectic eigenvalues are the two 
positive roots £1 > C2 > of the biquadratic equation 

C 4 -(det7 A + det7 B -2detcr Ai 3K 2 +det7 Ai 3 = 0. (15) 

The solution of Eq. (fT~5f) yields 

(I = \{D- v/£ 2 -4dct 7 AB) , (16) 



where D = det 7a + det 73 — 2&etoAB- The two-mode 
Gaussian state is entangled if and only if (2 < 1. In this 
case, we have Ej\f = — log(C2), while £V = otherwise. 
The condition (2 < 1 implies the necessary and sufficient 
entanglement condition D > 1 + det jab 0, , which 
explicitly reads, 

det 7a + det jb — 2 det oab > 1 + det 7,4s- (17) 

With the use of the method proposed in the present 
paper we can measure det Jab, det 7a, and det 73. An 
upper bound on det (tab in terms of these determinants 
can be derived from the condition that the symplectic 
eigenvalues Q of the covariance matrix jab must be 
greater or equal to one The lower eigenvalue £2 

is given by Eq. . where D is replaced with D' = 
det j a + det jb + 2 det oab- The condition Q > 1 yields 

2 det oab < det jab + 1 — det ja — det jb- (18) 

This, in turn, implies an upper bound on the lower sym- 
plectic eigenvalue £2 of the covariance matrix of . On 
inserting the upper bound on 2 det a ab given by Eq. I|18(l 
into Eqs. (|16fl and (|17fl we find that C2 < 1, so that the 
state is entangled when 

det j a + det jb > 1 + det jab (19) 

holds. Inequality (fH^I is thus a sufficient condition for 
entanglement, but it is not necessary as some Gaussian 
entangled states are not detected by this test. The main 
advantage of this test is that all determinants appearing 
in Eq. (|19[) can be determined by local measurements 
supplemented with classical communication between A 
and B. Moreover, if we can a priori assume that £ = 0, 
then measurements on a single copy of pab suffice. 

If we now want to exactly determine Ej\f, we also need 
a scheme to measure det gab ■ This can be accomplished 
provided that we can perform joint measurements on sev- 
eral copies of the state pab- Unlike the previous one, this 
scheme requires joint nonlocal measurements on modes 
A and B, and it involves several steps as schematically 
illustrated in Fig. 2. Using the scheme of Fig. 2(b), we 
can measure the determinants of the covariance matrices 
7 + and 7_ of modes A + and A_ that are linear combi- 
nations of the modes A and B, a± = (a ± b)/V2. After 
a simple algebra, we find that 

det(o- A B + o T ab) = 2 det 7+ + 2 det 7_ - det(j A + Jb)- 

(20) 

In order to determine det (ja + Jb), we have to carry 
out a joint measurement on two independent copies of 
the two-mode state, as depicted in Fig. 2(c). By mix- 
ing the modes A\ and B2 on a balanced beam splitter, 
we prepare an output single-mode state with covariance 
matrix ja+b = (ja + 1b) /2. Recall that the modes 
A\ and B 2 belong to two independent copies of the two- 
mode state pab, hence A\ and B2 are uncorrelated. Af- 
ter the measurement of det j + , det 7_ , and det Ja+ b , 
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FIG. 2: Direct measurement of the intermodal correlations 
det (tab of a two-mode Gaussian state, (a) Preparation of 
states with zero displacement £ = needed for steps (b) and 
(c). (b) Measurement of det 7+ (det 7_ is measured similarly), 
(c) Measurement of det7A+s- 



we calculate det {cab + &ab) f rom Eq. I|20() . It holds 
that det (cab + &ab) — 4 deters and the equality is 
achieved when a ab is symmetric. The matrix a ab can 
be brought to a symmetric form by applying a local phase 
shift exp(i(/>&^6) to the mode B using the phase shifter PS 
in Fig. 2. This transforms (tab to <jabU{4>) where 



U{(j>) = 



COS(p 

— sin< 



sin ( 
cos< 



To determine the phase shift <f> that symmetrizes (tab, 
we measure = Aet[a abU {(jr) + U T {4>)o~'ab\ f° r three 
different phases <f> = 0, 7r/4, and tt/2. It can be shown 
that 

V4> = Vo co s 2 <P + Vtt/2 sin 2 (j) + sin(20), 

where y„./4 = - (y + y^/ 2 )/2. The value of det a A B 
can be found as the maximum of yj, over (p which yields 



1 

det (tab = 77 



yo + j/tt/2 + J(yo- y^/2) 2 + 4y 2 /4 



This finally provides all the information required for the 
exact calculation of the logarithmic negativity Ej^. 

In summary, we have proposed a scheme for the direct 
measurement of the squeezing, purity, and entanglement 
of Gaussian states that does not require homodyne de- 
tection but only needs beam splitters and photodetectors 
with single-photon sensitivity. The scheme generally re- 
quires joint measurements on two copies of the state, but 
single-copy measurements suffice if it is a priori known 
that the mean (or coherent) values of the quadratures 
vanish, which is, e.g., the case of the squeezed and entan- 
gled states generated by spontaneous parametric down- 
conversion. We have shown that, based on Eq. lj*H3|) . the 
present method can be used to assess the entanglement 
of 1 x 1 Gaussian states by means of local measurements, 
without employing any local oscillator or interferometric 



schemes. Given the simplicity of the suggested setup, the 
prospects for an experimental realization in a near future 
look very good. 

Note added: The sufficient condition on entanglement 
[Eq. (|19Jl ] has recently and independently been derived by 
Adesso et al. 0] . Besides the lower bound on E^f linked 
to Eq. (|19J) . an upper bound on Etf expressed in terms 
of the determinants of 7^, 75, and 7,43 was also derived 
in Ref. |l8| . Remarkably, these two bounds are typically 
very close to each other, so the knowledge of the deter- 
minants of the covariance matrices provides quite precise 
quantitative information on the entanglement, making 
the direct measurement procedure particularly powerful. 
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